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Abstract 

Recently a quantization scheme for the phenomenological Maxwell theory of the 
full electromagnetic field in an inhomogeneous three-dimensional, dispersive, and 
absorbing dielectric medium has been developed and applied to a system con- 
sisting of two infinite half-spaces with a common planar interface (H.T. Dung, L. 
Knoll, and D.-G. Welsch, Phys. Rev. A 57, in press). Here we show that the 
scheme, which is based on the classical Green-tensor integral representation of 
the electromagnetic field, applies to any inhomogeneous medium. For this pur- 
pose we prove that the fundamental equal-time commutation relations of QED 
are preserved for an arbitrarily space-dependent, Kramers-Kronig consistent per- 
mittivity. Further, an extension of the quantization scheme to linear media with 
bounded regions of amplification is given, and the problem of anisotropic media 
is briefly addressed. 



1 Introduction 



Quantization of the phenomenological Maxwell theory of the full electromagnetic 
field in an inhomogeneous three-dimensional, dispersive, and absorbing dielectric 
medium of given permittivity necessarily requires a concept that is consistent with 
the principle of causality and the dissipation-fluctuation theorem, and necessar- 
ily yields the fundamental equal-time commutation relations of QED. Recently 
it has been shown [Q that the classical Green-tensor integral representation of 
the electromagnetic field in a medium with space-dependent, complex permittiv- 
ity can be quantized, in agreement with the conditions mentioned, introducing 
operator noise current and charge densities and expressing them in terms of a 
continuous set of bosonic fields. The quantization scheme generalizes previous 
work on dispersive and absorbing bulk material and one-dimensional slab-like 
systems with stepwise constant, complex permittivity 0, 0> S S S ■ 

In particular, the ordinary vacuum QED is recognized in the limit when the 
permittivity approaches unity, and the frequently used approximate quantization 
schemes for radiation in dispersionless and lossless inhomogeneous media (see, 
e.g., pi, |9], [TO]) and purely dispersive media (see, e.g., [|11], [12], |13|) are recognized 



in the narrow-bandwidth limit. Further, the concept is in full agreement with 
the Huttner-Barnett approach |TJ| to quantization of the electromagnetic field in 
bulk material. In this scheme, which is based on the Hopfield model of a ho- 
mogeneous dielectric, the electromagnetic field is coupled to a harmonic-oscillator 
polarization field that interacts with a continuous set of harmonic-oscillator reser- 
voir fields. All couplings are assumed to be bilinear, and the Hamiltonian of the 
total system is diagonalized by using a Fano-type technique [|1(| . 



The proof of the consistence with QED of the quantization scheme developed 
in |]]] requires the calculation of some frequency integral of the (classical) Green 
tensor in order to verify the fundamental equal-time commutation relation be- 
tween the electric and magnetic fields. So far the proof for a three-dimensional 
inhomogeneous medium has been based on the explicit expression of the Green 
tensor for a system that consists of two dispersive and absorbing bulk dielectrics 
with a common planar interface. Although it is the simplest inhomogeneous sys- 
tem, the involved form of the Green tensor needs performing a rather lengthy 
calculation, and the question about the validity of the theory for more com- 
plicated three-dimensional systems may arise. In this paper we show that the 
fundamental equal-time commutation relation between the electric and magnetic 
fields is satisfied for any inhomogeneous three-dimensional, dispersive and ab- 
sorbing dielectric medium, without making use of a particular form of the Green 
tensor. This enables us to show that the theory applies to arbitrary inhomo- 
geneous, linear media including media with bounded regions of amplification. 
Finally, we briefly address the extension of the theory to anisotropic media. 

The paper is organized as follows. The quantization scheme is outlined in 
Sec. 0. In Sec. R from the partial differential equation for the Green tensor 
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an integral equation is derived, and general properties of the Green tensor are 
studied. Section |] presents the proof of the fundamental commutation relation 
between the electric and magnetic fields, and in Sec. |^ it is shown that the scheme 
also applies to media with both absorption and (in bounded regions of space) 
amplification. Finally, a summary and some concluding remarks are given in 
Sec. |. 

2 Quantization scheme 

Following [|TJ, we spectrally decompose the (Schrodinger) electric and magnetic 
field operators as 

oo 

E(r) = J duE(r,uj) + H.c. (1) 

Q 

and 

oo 

B(r) = fduB{r,uj) +H.c, (2) 
o 

respectively, where E(r, u) and B(r, ou) satisfy Maxwell's equations 

V-B(r,cj) = 0, (3) 

V- [e e(r, w)E(r, u)] = p(r, u), (4) 
VxE(r,w) = «wB(r,w), (5) 

V x B(r,u) = -i—e(r,u)E(r,u) + /ioj(r,uj). (6) 

Here, the complex-valued permittivity 

e(r, u) = e R (r, u) + iei(r, u) (7) 

is a function of frequency and space, with 

e(r, a;)— >1 if to — > oo. (8) 

For chosen r the real part (responsible for dispersion) and the imaginary part 
(responsible for absorption) are related to each other according to the Kramers- 
Kronig relations, because of causality. This also implies that e(r,u) is a holo- 
morphic function in the upper complex frequency plane, 

d 

e(r,w) = (u;i>0). (9) 
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The dependence on r of e(r,io) indicates that the dielectric properties spatially 
change in general. 

In order to be consistent with the dissipation-fluctuation theorem, in Eqs. (^) 
and (§), respectively, an operator noise charge density p(r, lo) and an operator 

noise current density j(r, u) have been introduced, which fulfill the equation of 
continuity 

V • j(r,w) = itop(r,u). 



Eventually, j(r, u) is related to a bosonic vector field f(r, uj) as 



j(r, u) = lu \ — ei(r, lo) f (r, u) 



7T 



/i(r,o;), f](r',cu') = 5 l3 5(v - r')S(cu - u/), 







j?(r,u;),;j(rV; 



(10) 



(11) 

(12) 
(13) 



The fields f (r, u) for all u can be regarded as basic variables of an overall system 
that consists of the electromagnetic field, the polarization field and the reservoir 
fields and whose Hamiltonian reads 



(14) 



The quantization scheme implies that the electromagnetic field operators can 
be expressed in terms of f (r, to). From Maxwell's equations it is seen that E(r, to) 
satisfies the partial differential equation 



V x V x E(r,w) 



tO 



e(r, u;)E(r, to) = z/i ^j(r, u), 



(15) 



so that 



Ei(r,u) =ifx / d 3 suG ik (r, s,u)j As,u) 



(16) 



where j_ (s,u) is given by Eq. ([□]), and Gik(r, s, to) is the tensor-valued Green 
function of the classical problem. Here and in the following we adopt the conven- 
tion of summation over repeated vector-component indices. Combining Eqs. (|5|) 
and flTBj), the corresponding expression for Bi(r,u) is easily derived. The inte- 
gral representations of Ei(r) and Bi(r) are then found from Eqs. ([!]) and (Q), 
respectively, from which the (equal-time) commutation relations are derived to 
be 



EAr),Eu(r') 



5 t (r),5 fc (r') 







and 



Ei(r),B k (r')] = — e kmj d r ^ I du ^Gy(r, r', u) 



(17) 
(18) 
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(e/cmj, Levi-Civita tensor; d^^d/dx^). On the other hand, from QED it is well 
known that 

'E i (T),B k (r') 

which reveals that the quantization scheme is in full agreement with QED, if in 
Eq. fll8|) the integral over to yields 



+00 




—00 



It should be pointed out that this is also the condition for obtaining the correct 
commutation relations for the potentials and canonically conjugated momenta. 

Apart from scalar electrodynamics for slab-like systems |2|, |3], |], ||, |], [7], 
Eq. (|20D has been proved correct for bulk material and two infinite half-spaces 
with a common planar interface 0] by making use of the explicit form of the 
Green function. In what follows we show that the quantization scheme yields 
the correct commutation relations for arbitrary inhomogeneous dielectrics, i.e., 
for any permittivity e(r, u). For this purpose, let us first consider some general 
properties of the Green function. 



(19) 



3 Green function 



From Eqs. (|I5|) and flTop it is easily seen that the tensor-valued Green function 
[matrix elements of the fundamental solution of Eq. ([15])] satisfies the equation 



dldl 



5 ik A r + ^£(r,w) 



G kj (r,s,u) = 5ij5(r - s) 



(21) 



This equation and the boundary condition at infinity determine the Green func- 
tion uniquely. Similarly to Eq. (15) [together with Eq. fliTf)] there are no nontrivial 
solutions of the homogeneous problem. Let us consider absorbing bulk material. 
Since the Green function must vanish at infinity, absorption obviously prevents 
one from constructing a solution of the homogeneous equation which is different 
from zero at finite space points. When the dielectric material extends only over 
a finite region of space, we may assume that e(r, uS) — > 1 for 00. To preserve 
the analytical properties of e(r, u), the limit r — > 00 must be performed first, thus 
keeping a (small) imaginary part ei(r, u) in the permittivity, which again implies 
that there is only the trivial (zero) solution of the homogeneous equation. 

It is well known that the Fourier transform of a response function that de- 
scribes a causal relation between two physical quantities is a holomorphic function 
in the upper complex frequency half-plane (see, e.g., |T7|, [T8|, |T^]). A typical ex- 
ample is the causal relation between the averages of polarization and electric-field 
strength. Obviously, Z) i7 (r, s, cu) = ifj, u;Gij(r, s, to) as a function of to is nothing 
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but the Fourier transform of the tensor response function -Dy(r, s, r) that causally 
relates the electric field Ei(r,t) observed at space-point r and time t to an ex- 
ternal (point-like) current jj xt (r, t) = JJ xt (s, t)5(r — s) at space-point s and time 
t-r (r>0) [cf. Eq. ©]: 



£i(r,t) 



dTDij(r,a,T)Jj*(a,t-T). 



(22) 



Hence, 



i^ Q uGij{T,s,uS) 



f 

Jo 



(ire !Wr Ai(r,s,r) 



is a holomorphic function of to in the upper complex half-plane, i.e. 





duo* 



uG kj (r,s,uj) = (wi > 0), 



(23) 



(24) 



with 



wGfcj(r, s, to) — > if \u\ 



oo. 



(25) 

Note that Eq. ([24]) is in full agreement with the differential equation (|2T|) [together 
with the boundary condition at infinity]. In this equation the frequency uo is a 
parameter, and we may assume that Gy(r, s, o>) as a function of u; is differentiable 
with respect to uo in the upper complex half-plane. Applying 3/ duo* to Eq. (|2l"l) . 
we easily see that 
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dld r k -6 lk [A r + —e(r,uj 



_d_ 

duo' 



uoG kj (r,s,uo) = (cJr>0), (26) 



because of Eq. (§). Since there is no nontrivial solution of the homogeneous 
problem, we see that ooG k j(r, s, u) satisfies the Cauchy-Riemann equations 
From the theory of partial differential equations it is known (see, e.g., 
that there exists an equivalent formulation of the problem in terms of an integral 
equation. As shown in App. [A|, Gij(r,s,uo) satisfies the integral equation 



G i: j(r,s,u) = G^ 0) (r,s,u;) + / d 3 v K ik (r,v,uo)G kj (v,s,uo) 



(27) 



where 



and 



Gij(r,s,uj) = dij-d-d'q 2 (s,uo) g(\r-s\,u) 



K ik (r,v,u) 



^mg 2 (v,„) [^(|r-v|,o;)] 



+ 



g 2 (v,w) -ql{u) 



^(|r-v|,w) 5, 



ik- 



(2i 



(29) 
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Here, the function 



^ (|rU) = "i^r = J f^im (30) 



is introduced, where 



and 



q 2 (r 1 uj) = —e(r 1 uj) (31) 



c 2 



^ 2 



g » = ^e H, (32) 



c 



eo(c<j) = e(r, cu) being an appropriately space-averaged reference permittivity [for 
the integral equation with an s-dependent reference permittivity e (s, u), see 
App.0. Obviously, Gg } (r,s, uj) is the Green function for a homogeneous medium 
with permittivity e(r, u) = € (uj). The second term on the right-hand side in 
Eq. ([27]) essentially arises from the inhomogeneities. Note that according to 
the Fredholm alternative the solution of the integral equation (|27|) is unique, 
because of the non-existence of nontrivial solutions of the homogeneous problem. 
From Eqs. (|30D - ([32] ) and Eq. (§) it follows that the integral kernel Kik(r, v, lu), 
Eq. (|29|), is a holomorphic function of uj in the upper complex half-plane, with 

Kik(r, v, oj) — >• if \u\ — > oo, (33) 

where Kik(r, v, w) decreases as does e(r, cj) — 1. 

Let us write the integral equation (pTj) in the compact form 



where 



G = G (0) + ICG, (34) 



/CG = (/CG) y (r,s,u;) 

= y d 3 vK ik (r, v,uj)G kj (v,s,uj). (35) 

Assuming that G can be found by iteration, we may write 

oo 

G = G (0) + ]T/C n G (0) . (36) 



n=l 



From Eq. fl28|) it is seen that G^\r,s,u) has a cubic singularity |r — s| 3 for r- 



s, and Eq. (|29| ) reveals that the kernel ifjfc(r, v, a;) is only weakly singular (the 
singularity is weaker than the spatial dimension). Hence, at least after the third 
iteration step the result is perfectly regular at r = s. 
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4 Commutation relation 

The results given in Sec. ^| now enables us to prove Eq. (^0|) for arbitrary inho- 
mogeneous dielectrics. For this purpose we first decompose the Green function 
into two parts, 

Gij(r, s, u) = (Gi)ij(r, s, u) + (G 2 )ij(r, s, u), (37) 
where (Gi)^(r, s, uj) and (G2)ij(r, s, uS) satisfy the integral equations 

G lx = Gf+KG ll (^ = 1,2), (38) 

with [cf. Eq. (H)] 

{G 1 )f){r,s,u) = 8 ij g{\r-s\,Lo) (39) 

and 

(G a )g ) (r,s l w) = -d'dlq-^uMlr-slu). (40) 
It is easily seen that (G 2 )ij(r, s, u) can be given by 

(G 2 ) lJ (r,s,^) = ^r i (r,s,^), (41) 

where V is the solution of the integral equation 

r = r (0) + /cr, (42) 

with 

rf (r, s, u) = -dl q~ 2 (s, u)g{\T-a\, w). (43) 

Both ^(Gi)^^, s, u) and ^(Gy^ 1 *, s > a; ) are holomorphic functions of in the 
upper complex half-plane, with u(G ll )ij(r,s,u) — > if |a>| — > oo. Note that 
ui(G2)ij(r, s, a;) may be singular at u = 0. Nevertheless, when substituting Gij(r, s, ui) 
from Eq. ([37]) [together with Eq. (|4])] back into Eq. (|i~6D, we can integrate by 
parts and use the equation of continuity ( JT0| ) to obtain 

Ei(r,u) =i/x / rf 3 s^(G 1 ) iA .(r,s,^)j fe (s,^) 

+ /i y" d s sw 2 r j (r ) s ) w)£(B,w). (44) 

Hence, i[ioUj(Gi)ik(r, s, w) and /i ^ 2 rj(r, s, cu) are the Fourier transforms of the 
response functions relating the electric-field strength to the (noise) current density 
ju(s,u) and the (noise) charge density p(s,u) separately. Obviously, u 2 Ti(r,s,u!) 
is not singular at u = 0. 

Combining Eqs. (|37D and (|4l|), we easily see that the term on the right-hand 
side in Eq. flT8] ) can be rewritten as, on recalling that ekmjd^dj (. . .) = 0, 

+oo +oo 

ekmjd 7 ^ J du ^Gij(r,r\u) = e kmj d^ J du ^(Gi)y(r, r',o>)- (45) 
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Thus, only the noise-current response function ~ uj(G\)ij(r, r', uj) contributes 
to the commutator (|l"8"D . We now substitute in Eq. (|45|) for {G\)ij{v, r', uj) the 
integral equation (|38|) (/i = 1) to obtain 

+00 

/uj 
duj —(Gi)ij(r, r', uj) = iir5ij8(r - r ) 

—00 

+00 

+ J duj J d 3 v^K ik (r,v,uj)(G 1 ) kj (v,r',uj), (46) 

—00 

where the (bulk-material) relation |2| 

+00 

J duj £ (GOgV, r', oj) = iwS^ir - r') (47) 

—00 

has been used. Hence it remains to prove that the second term on the right-hand 
side in Eq. (|46|) vanishes [compare Eqs. (f45|) and (j^) with Eq. ([20])]. 

Since K ik (r,v,uj) and uu(Gi) k j(~v, s, uj) are holomorphic functions of uj in the 
upper complex half-plane, the uj integral can be calculated by contour integration 
along a large half-circle (with |u;| = R, R — >• 00). To calculate this integral, we 
recall that for |u;| ^00 both K ik {r, v, uj) and uj{Gi) k j{y, s, uj) approach zero at 
least as uj' 1 , and ujK ik (r, v, u;)(G ! i)fcj(v, s, uj) approaches zero at least as uj' 2 . 
Hence, for R—>oo the contour integral vanishes at least as -R -1 , and the second 
term on the right-hand side in Eq. fl4"6j ) indeed equals zero, i.e., 

+00 

/uj 
duj — {Gi)ij (r, r', uj) = iirS^^ - r), (48) 

—00 

which [together with Eq. (|45|) 1 shows, that Eq. (p0[ ) is valid for arbitrary space- 
dependent permittivity. In other words, the application of the quantization 
scheme to arbitrary inhomogeneous dielectrics yields the correct QED commuta- 
tion relation ([19]). 

5 Extensions of the quantization scheme 

The developed concept of quantization of the electromagnetic field in a disper- 
sive and absorbing background medium that is described in terms of a spatially 
varying, complex permittivity essentially rests on the following assumptions and 
principles, (i) The permittivity e(r, uj) of the dielectric medium and the Green 
tensor Gij(r, s, uj) of the classical Maxwell equations are holomorphic functions of 
uj in the upper complex frequency half-plane, because of causality, [ii) There is 
no nontrivial solution of the homogeneous Maxwell equations which satisfies the 
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boundary condition at infinity, i.e., the electric and magnetic fields are uniquely 
determined by their integral representations, (Hi) To be consistent with the 
dissipation-fluctuation theorem, noise current and noise charge densities must 
be introduced even if there are no additional sources embedded in the dielec- 
tric medium, (iv) Quantization then requires the integral representations to be 
regarded as relations between operator-valued fields, where the operator noise 
current density satisfies the commutation relation 



uj 2 — ei(r, u;)&,<J(r - r')5(uj - uj'). (49) 

7T 



So far we have assumed that ei(r, uj) > (uj > 0), as it is the case for absorbing 
media. Obviously, the statements (i) - (iv) remain valid, if (in agreement with the 
Kramers-Kronig relations) ei(r, uj) <0 (uj>0) in a bounded region of space, which 
corresponds to the presence of an amplifying medium in that region. Obviously, 
the commutation relation (|49|) is obtained if (in that region) the operator noise 
current density j(r, uj) is related to the bosonic field f(r, uj) as 



j(r,uj) = uj \ -— ej(r,uj)P(r,uj), (50) 

V 7T 

which reflects the well-known fact that amplification requires the roles of the 
noise creation and destruction operators to be exchanged (see, e.g., [^T], 
From inspection of Eqs. ([□]) and (|50|), we see that the two equations can be 



combined to express the operator noise current density associated with damping 
and amplification in terms of the bosonic field as 




i(r,uj) = ujJ^\e l (r,uj)\ 0( ei ) f (r, a;) + 6(-ei) ft(r, uj) , (51) 



with Q(x) being the unit step function [0(x) = l for x > 0, and Q(x)=0 elsewhere]. 
Note that the operator noise charge density p(r,u>) is given by Eq. (|T0|) , with 
j(r, uj) from Eq. (|5l|). 

The quantization scheme based on Eq. ( |5T| ) can be regarded as the extension of 
the concept for amplifying, one-dimensional slab-like systems p3| , |24|| to arbitrary 
inhomogeneous media that contain bounded regions in which amplification is 
realized. From the derivation given in Sec. ^ it is clearly seen that the fundamental 
QED commutation relation (|T^) is satisfied independently of the sign of ei(r, uj). 
For an absorbing medium the poles of u(Gi)ij(r, s, uj) as a function of uj are in 
the lower complex half-plane. When the gain owing to amplification (e.g., in 
a resonator-like equipment) tends to compensate for the losses, then the poles 
may approach the real axis and sharply peaked resonances may be observed. 
Obviously, if there are poles on the (real) axis, the uj integral must be performed 
along the axis uj + is, e — > 0. Note that in such a case the model of linear 
amplification may fail, because of nonlinear saturation. 
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Another possible extension of the quantization scheme is the inclusion in the 
theory of anisotropic media, for which the permittivity is a symmetric, complex 
tensor function of u, 

6y(r, uj) = eji(r,u), (52) 

which also varies with space in general. As we will show in a forthcoming article, 
the quantization scheme also applies to the electromagnetic field in anisotropic 
Kramers-Kronig dielectrics. The calculation relies on a symmetry relation satis- 



fied by the Green function according to the Lorentz reciprocity theorem [25, 26 



and the same integral relation ( |T8"D for the fundamental commutator between 
the electric and magnetic field operators can be derived. 



6 Conclusions 

We have studied quantization of the full electromagnetic field in linear, isotropic, 
inhomogeneous Kramers-Kronig dielectrics, using the formalism of Green-tensor 
integral representation of the electromagnetic field, in which the electromagnetic 
field operators are related to bosonic basic fields via the Green tensor of the classi- 
cal problem. The formalism can be regarded as the natural extension of the mode 
concept which only applies - apart from vacuum QED - to narrow-bandwidth 
fields. Basing on very general properties of the (classical) Green tensor, we have 
shown that the formalism yields exactly the QED equal-time commutation re- 
lations between the fundamental electromagnetic fields for any linear, isotropic 
dielectric medium. 

For this purpose we have derived an integral equation for the Green tensor, 
the kernel function of which describes the effect of spatially varying permittivity. 
From the holomorphic properties of the Green tensor and the integral kernel 
as functions of frequency it then follows that the QED equal-time commutation 
relation (|I9| ) between the electric and magnetic fields is preserved, independently 
of the dependence on space of the permittivity. Since the holomorphic properties 
are observed for absorbing media as well as amplifying media, the quantization 
scheme applies to any linear, isotropic, causal medium. The only condition is that 
amplification, which gives rise to a negative imaginary part of the permittivity, 
extends over bounded regions of space - a condition that is physically always 
fulfilled. It is worth noting that the scheme can also be extended to anisotropic 
media, as it will be shown in a forthcoming paper in detail. 

In order to show that the quantization scheme is consistent with QED, we 
have restricted attention to the equal-time commutation relations. Clearly, the 
results can also be used for determining the commutation relations of the (Heisen- 
berg) electromagnetic field operators at different times. Recalling that f(r, u,t) 
= f(r,u)e~ iwt , it can easily be derived that inclusion of cos[u;(t — t')] in the inte- 
gral on the right-hand side of Eq. ([18]) yields the commutator [Ei(r, t), Bk(r', t')]. 
Decomposing the Green tensor as shown in Eq. ( P?| ) [together with Eq. (fHP], we 
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find that 



j_ +00 

—00 

where (Gi)jj(r, r', a;) satisfies the integral equation (|38l), with (Gi) y (r, r', w) from 
Eq. (|39"D. This result can be regarded as the natural generalization of the well- 
known result of vacuum QED (see, e.g., |p8fl ). 
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= — tkmjd^ / duj — (Gi)y(r, r', u) cos[u(t-t% (53) 



A Derivation of the integral equation 



In order to derive an integral equation equivalent to the differential equation fl2T|) , 
we formally write in Eq. (|2l|) 

e(r, 10) = e(r, 10) + e Q (s, lo) - e (s, ct>), (54) 

where 6q(s, u;) is an appropriately chosen reference permittivity, which also satis- 
fies the Kramers-Kronig relations. Hence we may rewrite Eq. (]21f) as 



A r + q 2 (s,u) Gij(r,s,uj) 



Gij(r,8,u) 



+ d^d r k G kj (r,s,uj) - 5 ij 5(r - s), (55) 



where the abbreviations 



u 2 to 2 



have been used. Now we introduce the Green function 

e iq (s,w)\r\ 



9(\r\,s,u) 



47r|r| 



(56) 



(57) 



[q (s,u>) = (oj/c)Je (s,oj)\, which is easily proved to satisfy the differential equa- 
tion 

~A r + q 2 (s,iu)}g(\v\,s,uj) = -5(r). (58) 

The Green function gf(|r|,s, u) enables us to convert Eq. flSBP into the integral 
equation 

Gij(r,s,uj) = -J d 3 vg(\r-v\,s,u){ [q%(a,w) - <? 2 (v,u;)] G y (v, s,u) 

+ d?d v k G kj (-v, s, u) - 6^ - s)}. (59) 
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Next we apply d[ on Eq. fl2"l~| ) to obtain 



d-q (r,aj)G i:j (r,s,uj) = d^S(v-s), 

from which we find that 

dZG ij (r,s,uj) = -q- 2 (r,u;)d r j S(r- S ) 

9[lng 2 (r,cj) GL-(r,s,a>) 



(60) 



(61) 



Substituting in Eq. (|59f) for d^d k G k j(v, s, cu) the result of Eq. (|6l|), integrating 
by parts, and performing the integrals with (^-functions, we derive 



Gij(r,s,uj) = GP(r,s,uj) + J d\ K ik (r, v, s, u)G kj (v, s, w). 



(62) 



Here, 



Gf (r,s,^) = ^-9[%- 2 (s,u;) 



(ft r-s ,s,wj 



+ Si <9-g (v,u;)^(|r-v|,s,£j) 



(63) 



and the integral kernel reads 



K ik (r, v, s, uj) 



d v k hiq 2 (v,uj) [<9[#(|r-v|,s,^)] 



+ 



q 2 (v,u) - q 2 {s,uj) 



g(\r-v\,s,uj)5 ih 



(64) 



It should be pointed out that the reference permittivity e (s, uj) can be chosen 
freely in principle, since the exact solution of the integral equation (|B"2l ) does 
not depend on e (s, uj). In practice, however, it may be advantageous to choose 
e (s, uj) such that G^\r, s, uj) gives a sufficiently good zeroth-order approximation 
of Gij(r, s,uj) for an approximate solution of Eq. ([62]). 

In the simplest case eo(s,u) may be chosen to be independent of s, e.g., by 
averaging e(r, uj) over space, 



e (s,uj) -> e (uj) = e(r,uj) . 



(65) 



Obviously, in this case the Green funcion (|5T| ) and the integral kernel (^) become 
independent of s and Eqs. (|62]) - fl6~4]) reduce to Eqs. ( ^7|) - (|9|) [<?(|r|, s, a;) — > 
g(\r\,uj), K ik (r,v,s,uj)^>K ik (r,v,uj)] together with Eqs. ([ft]) - (|52|). 
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